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ABSTRACT 

Research questions that test the interaction between 
either categorical variables or multiplicative expressions are 
commonly found in multiple regression and analysis of variance 
applications. In contract, research questions posing an interaction 
between observed variables in a path analytic model are not commonly 
found in the literature. This is due in part to some of the issues 
and problems that arise in relational path models. Path models with 
latent variables also pose problems for the researcher that make the 
inclusion of interaction tests difficult. This paper further explores 
the categorical, group, or multi-sample approach and the continuous 
or multiplicative approach to testing interaction effects in 
structural equation models that use latent variables. Some software 
programs have provided a way to test simple interaction effects among 
latent variables. Examples are included. The developers of structural 
equation modeling software programs are encouraged to develop 
procedures for testing interaction effects. Until then, the testing 
of interaction effects in structural equation models will not be easy 
to perform. Of the two approaches discussed, multi-sample and 
multiplicative, the multi-sample approach is generally an easier 
technique to use for testing the equality of coefficients. An 
appendix contains a multi-sample interaction example and a LISREL 8 
program example. An additional attachment contains other software 
examoles and the outputs of some programs. One table and one figure 
illustrate the discussion. (Contains 29 references,) (Author/SLD) 
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ABSTRACT 



Research questions that test the interaction between either categorical variables or 
multiplicative expressions are commonly found in multiple regression and analysis of variance 
applications. In contrast, research questions posing an interaction between observed variables in 
a path analytic model are not commonly found in the literature. This is due in part to some of the 
issues and problems that arise in relational path models. Path models with late»:t variables also 
pose problems for the researcher that make ihe inclusion of interaction tests difficult. This paper 
further exolores the categorical, group or multi-sample approach and the continuous or 
multiplicative approach to testing interaction effects in structural equation models that use latent 
variables. Some software programs have provided a way in which to test simple interaction 
effects among latent variables, and examples are included in the paper. The developers of 
structural equation modeling software programs are encouraged to further develop procedures for 
testing interaction effects. Until then, the testing of interaction effects in structural equation 
model , will not be easy to perform. Of the two approaches, multi-sample and multiplicative, the 
multi-sample approach is generally an easier technique to use for testing the equality of 
coefficients. 
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TESTING INTERACTION EFFECTS IN STRUCTURAL EQUATION MODELS 



Research questions that test the interaction between either categorical variables or 
multiplicative expressions are commonly found in multiple regression and analysis of variance 
applications. In fact, the rationale for a 2 x 2 fixed-effects factorial study is to first test whether 
the interaction effect between variables is significant. In the presence of a significant interaction 
effect, the main effects are not interpreted. The basic reasoning behind testing for interaction is 
that two variables when interacting produce an effect different from their main effects. For 
example, hydrogen and oxygen as separate elements (main effects) have specific gaseous states, 
however, when they interact water is produced. Their interaction produces a different effect than 
their individual main effects. A ratiier simplistic example, but it serves to illustrate the basic 
nature of testing for interaction effects. Since researchers encounter both categorical variables 
and continuous variables, different approaches have emerged for testing interaction among the 
variable types. In this paper, tiie term multi*sample will be used when referring to categorical 
variables and the term multiplicative will be used when referring to continuous variables in 
testing interaction effects. 

Types of Interactions 

Interaction effects may take many forms. The interaction between two variables may 
enhance tiie effect which each variable has separately, as when alcohol and barbiturates are 
ingested together. In contrast, the interaction may dampen the individual effects of the two 
variables, as when two noises combine to create a zone of apparent quiet. Smith and Sasaki 
(1979) describe a 'Consistency'' interaction effect: when the two variables both deviate from their 
means in the same direction, their joint effect may be enhanced, but if the deviations are in 
opposite directions, the interaction effect will be opposite, as well. Also, two variables can 
interact forming a third variable which is entirely different from each variable's individual effect. 
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Multiple Regression and Analysis of Variance 

The testing of interaction effects in multiple regression and analysis of variance have 
become quite common, and covered in many earlier textbooks (Draper & Smith, 1966; McNeil, 
Kelly, & McNeil, 1975; Cohen & Cohen, 1975). The reader is probably aware of the various 
terminology used, namely, polynomial, curvilinear, or non-linear effects used to express a model 
with an interaction term. A model using continuous variables would take the form: 



where X, represents a multiplicative interaction, X,*X2 , or could include interaction effects 
expressed in a model as a polynomial term of the form: 



where X^i represents the square of the scores on X, More recent treatment of the topic by 
A:ken and West (1993) further illustrates examples using continuous predictors, categorical 
predictors, and interactions between categorical and continuous variables. 

In general, an interaction effect exists whenever the character of the relationship between 
variables A and B is affected by the level of variable C (Keppel & Zedeck, 19X9). We can 
demonstrate this by examining a multiple regression equation without interactions: 

A = p, *B + p2*^^ (') 
The first partial derivatives of this equation with respect to B and C are. respectively, P | and p^. 
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These derivatives show, for example, that the only effect of a change in B on variable A is the 
immediate effect through Pj. By contrast, consider a regression equation with an interaction 
term: 

A = pi *B + p2*C + P3*(B*C) (2) 

Now the first partial derivatives are pj + P3 * C and P2 + P3 * B, respectively. In other words, 
changes in C, for example, affect A not only directly but also indirectly by changing the character 
of the relationship between A and B. Accordingly, the estimated direct effect of B on A, p j , can 
be thought of as an "average" value computed over the range of values for variable C which are 
represented in the data set (Keppel & Zedeck, 1989). 

Path Analysis Models 

Path analysis provides an approach to testing models that indicate direct and indirect 
effects among specified variable relationships (Pedhazur, 1982; Pedhazur & Schmelkin, 1991). 
Certain assumptions are made when conducting path analyses beyond those inherent in 
conducting multiple regression analys(*s, and It is recommended that these be checked prior to 
using the technique (Newman, 1991 ). Other strengths and problems in conducting path analyses 
have also been elaborated by Pohlmann (1991). Research questions posing an interaction 
between observed variables in a path analytic model, however, were not commonly found in the 
research literature (Newman, Marchant, & Ridenour, 1993). Of the studies reporting a test of 
interai tion between variables, the most common approach to test interaction effects was to 
analyze separate models and compare model weights. In a few studies, a first-order interaction 
term was included. The authors concluded that, although interaction was implied or explicitly 
stated in the theory, the models in many of the studies failed to consider interaction effects. 
Pohlmann (1993) has further indicated that nonlinear relationships and latent variable assumptions 
can lead to serious specification errors in structural models. Path models with latent variables 
obviously pose problems for the researcher that make the inclusion of interaction tests difficult. 



INTERACTION IN LATENT VARIABLE MODELS 



Interaction effects in multiple regression and path analysis with observed variables has 
posed problems for researchers. These problems with corresponding assumptions are exacerbated 
in SEM. First, there is the specification problem. Linear models ease the task of determining 
which relationships to investigate and simplify distributional assumptions. Discarding the linearity 
restriction magnifies the critical role theory plays in focusing the structural equation modeling 
research effort. Researchers modeling interaction effects must also be sure to collect data which 
includes a range of values where interaction effects occur. Nonlinear functions may appear 
linear within certain ranges, so it seems prudent to suggest that observed variable values be 
plotted and investigated prior to model inclusion. In addition, the predictors might be normally 
distributed, but their joint-distribution not normally distributed in a multiplicative interaction 
term. The test of interaction effects are further complicated by the degrees of freedom 
introduced, subsequent interpretation of various fit indices and estimated standard errors, 
which are only approximate under certain conditions. In the multiplicative model when 
multiplying two observed variables, the product will contain the "true scores"and error with the 
distribution of the products of the error variables not normally distributed even if the error 
variables themselves are These concerns, as well as, several others discussed in the paper 
complicate testing interaction effects in structural equation models. 
There are two general approaches for dealing with interaction effects in structural equation 
models. Interaction ; involving discrete or discretized (categorical) variables can be evaluated 
using multi-sample analysis, or researchers can model interactions between two continuous 
variables by creating measures of the interaction ''construct," using techniques similar to those 
pioneered by Kenny and Judd (1984). Each of these techniques has its own strengths and 
weaknesses, and will be discussed separately. The multi-sample approach is more readily 
adaptable to testing the equality of factor structures, equal regression coefficients, and means of 
latent variables (Joreskog & Sorbom, 1^^93a). 
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Multi-Sample Approach 

In the multi-sample approach, the different samples are defined by the different levels of 
one or both of the interacting variables. If interaction effects are present, then certain parameters 
should have different values in the different samples. For example, in the context of an 
experiment, Mackenzie and Spreng (1992) tested an interaction model of advertising 
effectiveness. Their model proposed, in part, that a person's attitude toward an advertisement 
(Aj^) interacted with the person's motivation to process the ad to affect the person's attitude 
toward the advertised brand (A^). The structural model may be represented as: 

h\ = a, + Yii + Ci (3) 

where r| j is A^^, a j is the intercept, ^ j is Aj^, and Cj is an error term. Mackenzie and Spreng 
hypothesized that subjects with a higher motivation to process the ad would have a hi^iher mean 
level of after viewing the ad. This would be reflected in a higher intercept term for the 
subjects in the high motivation condition. They proposed that A^^ should have a positive "main" 
effect (as represented by y 1 1 ) on A^^. However, they further proposed that the interaction 
between A^^ and motivation would make the A^^-- Aj^ relationship ^i^akfil in the high motivation 
condition (see Figure 1). 



insert Figure 1 about here 
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Under the multi-sample approach, researchers can investigate interaction effects using x 

difference tests. In this example, the two samples are composed of subjects who were either 

exposed or not exposed to a manipulation designed to increase motivation. To test for the 

interaction between motivation and A^, first estimate a model where Y 1 1 restricted to be equal 

across the two groups, and then estimate a model where the parameter is allowed to differ in the 

2 

two samples. The two models are nested, so that, given assumptions, the test statistic is a x 
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difference test with one degree of freedom. In the example drawn from Mackenzie and Spreng 
2 

(1992), the X statistics for the two models were 44.55 (25 df) and 36.38 (24 df). The difference, 
8.17 (1 df) is significant (p < .005), so the interaction hypothesis is supported. The maximum 
likelihood estimate for y 1 1 was 0.881 (t = 13.99) for the low motivation subjects but only 0.644 
(t = 9.58) for the high motivation subjects. The Appendix contains the LISREL 8 command file 
for this problem. 

The multi-sample approach is the natural choice when the interaction relationship involves 
one or more categorical variables, and this approach has additional strengths as well. It can be 
used to represent a very wide variety of interaction effects without requiring substantial new 
methodological developments. Even the extension to higher order interaction effects is fairiy 
obvious. This approach can be used regardless of whether the interaction intensifies or mutes the 
effects of the individual variables. Because the interaction effect is represented in the difference 
between samples, the researcher may be able to preserve linear relations between variables within 
each sample, thus avoiding potentially significant complications in fitting the rr.cdel. Finally, 
nearly every commercially available SEM software package allows for multi-sample analysis with 
restrictions across samples, so researchers can probably use the package with which they are most 
familiar. 

The multi-sample approach also has a number of weaknesses. First, this method forces the 

researcher to divide the total available sample size by the number of groups. In a methodology 

which relies so heavily on asymptotic properties, and where those properties may only be 

achieved at relatively high sample sizes, this is a serious limitation. Furthermore, this reduction in 

2 

sample size may confound results from x difference tests. MacCallum, Roznowski and 

Necowitz (1992) noted a substantial degree of instability of fit indices in even moderately sized 

samples. Thus, it is possible that the multi-sample approach may yield subsamples which are so 

2 

small that instability '\ the X statistic will mislead the researcher into believing that an interaction 
effect exists, whether it does or not. On the other hand, the researcher may be able to minimize 
this problem by estimating only those distinct parameters that are necessary. In particular, 
researchers usually have no reason to expect that measurement parameters, such as loadings, will 
be different in the various subsamples. Specifying equality for such parameters may partially 
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lessen the impact of sample division on results. 

Additional problems arise if researchers attempt to apply the multi-sample approach to 
interactions involving two continuous variables. In such cases, the researcher must reduce at least 
one of the continuous variables to a categorical variable, for purposes of defining the groups. The 
categorized variable will contain less information than the original continuous variable, and loss of 
infi nation is undesirable (Russell & Bobko, 1992). Furthermore, in the process of 
categorization, a number of observations will probably be misclassified due to measurement error. 
SEM researchers account for measurement error by using multiple measures and latent variables, 
but latent variables are typically modeled as continuous. Perhaps researchers could compute 
factor scores and perform the transformation on that basis. This approach however has not been 
reported in the literature. Furthermore, choosing the point(s) at which to divide the sample into 
the various groups is(are) arbitrary. If researchers use arbitrary values (such as the median), 
mere random sampling error will ensure that many cases are misclassified, violating the basic 
assumption that the cases in a particular subsample all come from the same population. 

The Continuous- Variable Approach 

When an interaction involves two continuous variables, researchers may wish to adapt 
procedures similar to those described by Kenny and Judd (1984). In the continuous-variable 
approach, an interaction construct, for example ^3, is represented explicitly, along with ^ j and 
the "main effect" c(^nstructs, in a structural model of the form: 

ni = Yi *5i + 72*^2 ^ Y3*53 ^ ^1 ^"^^ 

where r\ | is the dependent construct and Ci is an error term. Notice that the relationship between 
r| I and ^3 is itself linear. Thus the researcher must create ^3 so as to represent the kind of 
nonlinear relationship believed to exist between T] | and ^ | and ^2- For example, to check for a 
multiplicative interaction between ^ I and ^3 '^^^t be made equal to * ^2' Measures or 
indicators of $3 are then created as functions of the measures of 5| and in a similar way. 

The structure of the interaction model emerges as a logical extension of the measurement 
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model for i\ and 1^ should be noted that these developments are due to the pivotal work by 
Kenny and Judd (1984). The basic facte r analysis model Cuii be stated as: 



(5) 



where x is a vector of manifest variables, A is a matrix of loadings, and 8 is a vector of 
measurement error terms. In thinking about the continuous-variable approach, it will help to keep 
' number of points in mind. Fir ;t, remember that both ^ and 6 are latent variables affecting x. 
While researchers may habitually draw sharp mental distinctions between the two, there is little 
difference between them, beyond the basic fact that one set of latent variables are common to 
several manifest variables while the other set load uniquely on individual manifest variables. 

Second, remember the assumptions underlying the basic factor model. The common 
factors, ^, are assumed to be uncorrected with the unique factors, 8. Both the common and 
unique factors are assumed to be normally distributed and have zero means. We can represent the 
covariance matrices of the common and unique factors as O and 6, respectively. In modeling a 
multiplicative interaction, Kenny and Judd (19S4) cited a third useful result from the theory of 
normaUy distributed variables. If A and B are normally distributed, then: 

2 2 2 2 

o (A*B) = o (A)*a (B) + (o(A,B)) (6) 

That is, the variance of the product is equal to the product of the variances plus the square of the 
covariance. Obviously, if A and B are uncorrected, then the last term disappears. 

From this point, Kenny and Judd (19S4) used simple algebraic substitution to develop 
their model of multiplicative interaction effects. Suppose ihcre are two measures which foUow 
factor models: 




(7) 



Then the product x-^ = Xj * X2 should be reflected in the following model as: 



X3 = X,*X2* 51*^2 + ^1*^1**2 ' ^2*^2**1 + *1*^2 



(9) 



Or 



X3 = X3 * ^3 




(10) 



where 




(11) 



(12) 



(13) 



(14) 



(15) 



All 3f thiese new latent variables are mutually uncorrelated and uncorrelated with all other latent 
variables in the model. 

In order to incorporate this interaction effect into a model, Kenny and Judd (1984) needed to 



relationships. This involved specifying some model parameters as nonlinear functions of other 
parameters. Kenny and Judd (1984) used a version of the COS AN program (Fraser, 1980) to 
accomplish this. Hayduk 0987), applying Rindskopfs (1984) work with "phantom" constructs, 
showed how the Kenny-Judd model could be specified artificially using the simple equality 
constraints found in programs like Joreskog and Sorbom's LISREL software. However, the 
difficulty of using Hayduk's method was sufficient to dissuade most researchers from attempting 
to specify such an interaction effect in a model. Today, both SAS/Proc CALIS (19^0) (cited by 
vValler, 1993; program by Wood in Appendix) and Joreskog and Sorbom's (1993b) LISREL 8, as 
well as, Amos and Mx software programs (see Appendix), include the ability to specify nonlinear 
constraints directly. 



specify X3 as a function of latent variables whose variances and covariances reflected these 
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LISREL8 EXAMPLE 
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With LISREL 8, researchers can specify the Kenny and Judd (1984) interaction model 
using V Alue and EQuality commands and the new Constraint commands. For example, the 
i .ny-Judd model implies that: 

o^(^3) = o^(5i)*o^(52) + (a(^,,^2>)^ (16) 
Researchers can specify this relationship using the Constraint line: 

CO PH(3,3) = PH(1,1) * PH(2,2) + PH(2,1 )**2 
Similarly, the model implies that: 

o^(^4) = o^(5,) * 0^(62) (17) 
and this relationship can be specified via: 

CO PH(4,4) = PH(M) * TD(2,2) 

One peculiarity of LISREL 8 is that all terms on the right side of a CO statement must be 
free parameters-they ct'.nnot be fixed values and they cannot themselves be subject to constraints 
(Joreskog & Sorbom, 1993b). If a parameter is fixed, then the researcher should replace the 
parameter in the Constraint line with the fixed value. If the parameter is constrained, then 
substitution may eliminate the problem. Furthermore, LISREL 8's routine for computing starting 
values will not work for the continuous-variable interaction model. The starting values chosen by 
the researcher can have a tremendous impact on the ability of the software to converge on a 
solution, even when the model is correct in the population. 

Kenny and Judd (1984) demonstrated their approach by fitting a covariance matrix 
computed from 500 Monte Carlo observations drawn from a population defined by an interaction 
model with known parameter values. A complete LISREL 8 program specification for the Kenny 
and Judd model is in the Appendix. Table 1 reports the true parameter values and the estimates 
produced by Kenny and Judd (1984), Hayduk (1987), and the LISREL 8 program in the 
Appendix. All three sets of results are very similar to the true parameter estimates, and to each 
other. 
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Insert Table 1 about here 



Researchers should consider using the continuous-variable approach particularly when 
wanting to avoid the weaknesses of the multi-sample approach which included: splitting the 
sample and possibly misclassifying observations; categorizing a variable and losing information; 
and increased sample size requirements. Kenny and Judd (1984) noted a possible benefit of using 
the continuous-variable approach, namely, parsimony. All but one of the additional parameters 
involved in the interaction model are exact functions of the "main effects" parameters. The only 
new parameter is the y parameter linking the interaction construct to the dependent construct. 

However, the continuous-variable approach has its own weaknesses, and specifying a 
continuous-variable interaction model is still a tedious, complex undertaking, with a high risk of 
making programming errors. Second, if a model includes several measures of a construct, or 
higher-oraer interaction terms, the technique becomes even more difficult to program. If 
constructs A and B have a and b measures, respectively, then the interaction construct AB will 
have a * b measures. If each construct has four measures, as is necessary to apply Anderson and 
Gerbing's (1988) internal consistency criterion, then the interaction construct will have 16 
measures. With the four measures of the two "main effect" constructs, and at least one measure 
of the dependent construct, the model will have at least 25 manifest variables before considering 
any other constructs. 

Under the continuous-variable approach, the researcher must specify the functional form 
for the interaction, construct measures of tlie interaction construct, and then model the 
relationships between the interaction measures and the latent variables. The Kenny-Judd model 
involved a specific type of multiplicative interaction, but this is hardly the only form that an 
interaction can take. For other types of interaction there is littie prior work available to guide 
researchers. Researchers may also face a profound niulticollinearity problem. It is very likely 
that the interaction measures will be highly correlated with the measures used to construct them. 
This niulticollinearity can be devastating to factor analytic measurement models, causing measures 
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to be more highly correlated with measures of other constructs than they are with measures of the 
same construct. For multiplicative interactions betw en normally distributed variables. Smith and 
Sasaki (1979) demonstrated that researchers can eliminate this multicoUinearity by centering the 
original variables-expressing them as deviations from their means-before computing the product 
variables. However, centering the variables alter.> .le form of the interactive relationship. Notice 
that the product variable A * B will have the same value when A and B are both well above their 
means and when they are both well below their means. Thus, Smitli and Sasaki (1979) noted that 
centering the variables turns the original multiplicative interaction into a kind of "consistency 
effect." Researchers who want to model other types of interactions may find no easy answer to 
the problem of multicoUinearity. 

The interaction model also presents distributional problems which are much more serious 
than those associated with SEM techniques in general. First, the formula for the variance of a 
product variable (Equation 6) assumes that the two original variables are normally distributed. If 
the original variables are severely nonnormal, the variance of the product variable can be very 
different from the value implied by Equation 5, and the interaction model will peiform poorly. Of 
course, suitable data transformations may result in approximate normal distributions for the 
original variables. There is, however, a second distributional problem. Many, if not most, 
functions of two normally distributed variables will not themselves be normally distributed. This 
problem applies to both the measures and the latent variables. This nonnormality violates 
distributional assumptions associated with estimation methods such as maximum likelihood. 
Furthermore, estimation methods that do not make distributional assumptions may not work for 
interaction models. Joreskog and Yang (1995) indicated that the asymptotic weight matrix 
associated with the covariance matrix for an interaction model may be not positive definite, due to 
dependencies between moments of different variables which are implied by the interaction model. 
Distribution-free estimation methods that use the inverse of the asymptotic covariance matrix will 
therefore fail. Researchers may minimize this problem by estimating the model using a 
distribution-sensitive method but evaluating the results with a bootstrapping technique. The 
rescaled bootstrapping method outlined by Bollen and Stine (1993) may be particulariy 
appropriate in this case. Note that techniques for computing bootstrap estimates of the weight 

lb 



matrix will not be helpful in this case, because the weight matrix is singular in the population 
(Yung & Rentier, 1994). 

CONCLUSION 

Clearly, both the multi-sample and continuous variable methods have their weaknesses. 
While new developments in SEM software are creating exciting new possibilities, researchers 
need to remain focused on their research hypotheses and theory. The multi-sample approach 
appears to be the more generally useful of the two techniques and gives rise to fewer problems, if 
theory, hypotheses, and variables are appropriately considered. Researchers who hypothesize 
"consistency effect" interactions are probably advised to apply the Kenny and Judd (1984) 
approach, but for other researchers, the difficulties of the continuous variable approach may not 
be worth the effort. Researchers may balk at the sample-splitting that is a part of the multi-sample 
approach, but in our experience, this problem is only a major concern in cases where researchers 
are working with data sets that were not collected with an interaction model in mind. Researchers 
will always face significant problems when they attempt to use data for unintended purposes, 
especially in SEM. Furthermore, it is not clear, in general, which of the two approaches described 
here is more robust to this kind of ex post facto modeling. Overall, we encourage further 
development of the software programs to more readily incorporate tests of interaction effects and 
address many of the issues and concerns indicated by the research literature. 
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MULTI-SAMPLE INTERACTION EXAMPLE 

(Annotated) 

LISREL8 conunand lines are flush left Comments are indented. 

When estimating multi-sample problems in LISREL, it is important to include a different title card 
for each group, m order to make the output more readable. 

TITLE CARD: Low Motivation Sample 

The DA card includes the NG parameter, which specifies the number of groups or 
samples. As in all multi-sampie problems^ researcTiers should be careful to analyze the 
covariance matrix rather than the correlation matrix (Joreskog & Sorbom, 19X9). 

DA Nl=6 NO=200 MA=CM NG=2 

This problem requires both covariances and means. These keywords would be followed 
by eiuier the moments themselves or by their file locations. 

CM 
ME 

LA; ATT_AD1 ATT_AD2 ATT_AD3P ANDATl BRANDAT2 BRANDAT3 

Mackenzie and Spreng (1992) specified their model in terms of t) constructs and y 
variables only. While it is not required for current purposes, it does simplify the 
programming. 

MO NE=2 NY=6 LY=FU.F1 BE=FU.FI PS=SY TE=SY TY=FR AL=F1 

LE; ATT_AD BRANDATT 

FR LY(2.1) LY(3.1) LY(5.2) LY(6,2) 

VA 1 LY(1,1)LY(4.2) 

FR BE(2,1) 

OU AD=OFF 

As Joreskog and Sorbom (1989) note, setting up multi-sample problems with LISREL is a 
less onerous task than one might expect. Specifications for the first group in the set-up 
become default values for later groups. 

TITLE CARD: High Motivation Sample 

DA NO=16() 

CM 

ME 

LA; ATr_ADl ATT_AD2 ATT_AD3 BRANDATl BRANDAT2 BRANDAT3 

The construct intercepts (the a's) are only defined relatively across the two groups. Thus, 
the intercepts for the first group are fixed to 0. Then estimates of the intercepts in the 
second group are evaluated relative to 0. 

MO LY=IN BE=PS TE=PS PS=rS TY=IN AL=FR 
OU 
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LISREL 8 PROGRAM FOR THE KENNY-JUDD MODEL 

Program lines are flush left. 
Comments are indented one-half inch. 

Title card: Kenny-Judd (1984) simple example with fewest starting values 

DA NI=9 NO=500 MA=CM 

Kenny and Judd labeled their independent constructs X and Z, and labeled their measures 
similarly. The underlines in these labels indicate multiplication. 

LA; xl x2 zl z2 xl_zl xl_z2 x2_zl x2_z2 y 

SE; 912345678 

CM 

For most parameter matrices, most elements will be fixed, so it makes sense to begin with 
these matrices set to Fixed. 

MO NY=1 NE=1 NX=8 NK=7 PS=SY TE=F1 PH=SY,FI TD=SY,F1 GA=F1 

This program includes ETA and KSl in the latent variable labels in order tc minimize 
confusion between measures and constructs. Note that, for the KSl latent variables, 1 and 
2 are the main effects, X and Z, while 3 is the interaction, X*Z. 

LE; ETA_Y 

LK; KS1_X KSI_Z KSLXZ X^TD3 X_TD4 Z_TD1 Z_TD2 

These fixed loadings set the scales of constructs Y, X and Z. 
VA 1.0 LY(K1)LX(1,1)LX(3,2) 

These are the free parameters relating to the main effect constructs and their measures. 

FR GA( 1 , 1 ) G A( 1 ,2) PH( 1 , 1 PH(2,2) PH(2, 1 ) 

FR LX(2,1) LX(4,2) TD(1,1) TD(2,2) TD(3,3) TD(4,4) 

The direct effect of the interaction construct on the dependent consruct is the only 
additional free parameter in the interaction model. 

FRGA(1,3) 

Following this are the constrained elements of the interact model. The loadings of the 
interaction measures on the interaction construct are simple products of the loadings of xl 
and x2 on X and zl and z2 on Z. Since two of those loadings are fixed to 1, they ao not 
appear in the Constraint statements. Remember that the parameters that actually appear 
on the right side of a Constraint line must themselves be entirely FRee. 

VA l.()LX(5,3) 

CO LX(6,3) = LX(4,2) 

CO LX(7,3) = LX{2,1) 

CO LX(8,3) = LX(2,1 ) * LX(4,2) 
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Note that this program puts the measurement error variances into 0t. Alternately, the 
program cojld have created more 5 constructs to represent these-bcft that would 
needlessly increase the dimension of several parameter matrices, slowing down the 
software. With 8 exogenous measures, 65 will be 8 x 8 anyway, so why not put the 
product measurement error variances there? Note that the covanances of these 
measurement error terms are all zero. That is why those covariances do not appear in 
these constraint lines. 

CO TD(5,5) = TDd.l ) * TD(3.3) 

CO TD(6,6) = TD(1,1) * TD(4.4) 

CO TD(7,7) = TD(2.2) * TD(3,3) 

CO TD(8,8) = TD(2.2) * TD(4.4) 

On the other hand, this program specifies that the constructs X and Z are correlated, so 
we must remember their covariance in specifying the variance of X*Z. 

CO PH(3.3) = PH(l.l) * PH(2,2) + PH(2.1)**2 

And these constraints represent the loadings of the interaction measures on the nuisance 
constructs, the products of constructs and error terms. In looking at these constraints, 
remember that, in LISREL 8, the loading of a measure on its own measurement error is 
implicitly fixed to 1. Also remember that LXd^l) and LX(3.2) are fixed to 1. All of those 
fixed I s are 'invisible' elements of these constraints. 

VA 1.0 LX(5,4) 

CO LX(7.4) = LX(2,1) 

CO PH(4,4) = PH( 1 , 1 ) * TD( 3,3 ) 

VA 1.0 LX(6,5) 

COLX(8,5) = LX(2,l) 

CO PH(5,5) = PH(l,r) * TD(4,4) 

VA l.()LX(5,6) 

COLX(6,6) = LX(4,2) 

CO PH(6,6) = PH(2.2) * TD(l.l) 

VA 1.0 LX(7.7) 

CO LX(8.7) = LX(4,2) 

CO PH(7.7) = PH(2.2) * TD(2,2) 

The Kenny-Judd model will run in LISREL 8 with only these starting values specified. In 
some cases, the problem will require starting values for virtually every parameter in the 
model. Do not ignore the constrained parameters, either. Just because the 'original' 
parameters have starting values. LISREL 8 will not automatically compute starting values 
for constrained parameters on that basis. The NS parameter on the output line suppresses 
automatic starting value calculation--this may eliminate the occasional processing error. 

.ST-.15GA(1,1) 

ST.35GA(1,2) 

.ST.70GA(1.3) 

LISREL 8's admissibility check will be unhelpful for interaction models, and shoukl be 
turned off. ML is Ll.iREL's default fit function, but it is sp -cified here for clarity. 

OU ML AD=OFF 

2''\ 
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Table 1 

Comparison of Parameter Rrimatfts 

Estimates 

Parameter Population Value Kenny-Judd (1984) Hayduk(1987) LISREL8 



Yii -.15 -.17 -.17 -.17 

Y21 -35 .32 .32 .32 

Y31 .70 ,70 .71 .71 

.60 65 .64 .65 

X42 .70 . 69 .69 .69 



f 



22 



Figure Captions 



Figure 1 . An example of an interaction effect. 
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Attitude toward 
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Low Motivation 




High Motivation 



Attitude toward 
the Ad (Aa) 



Yi : Main effect of Aa 

a2 - tti : Main effect of Motivation 

Y2 " Yi • Interaction effect 
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MX SOFTWARE EXAMPLES 



mod. sum: 

Moderator variables: two approaches. 

The standard approach allows for main effects of a predictor variable and a moderator 
variable, along with a multiplicative interaction term predictor*moderator. This may be modelled 
with either summary statistics such as means and covariances, or with the raw data. 

The Mx program, moda.mx, illustrates this procedure with simulated data on 500 
subjects. A structural equation model of multiple regression of the outcome (Y) variable on the 
three X variables is fitted by maximum likelihood to the raw data vectors. The last two lines of 
the script revise it to drop out the regression on predictor*moderator. The difference in the fit 
function between these two runs is 120.6 which is distributed as chi-square with 1 d.f. ( i.e., 
5270.752 - 5150.144= 120.6). 

The Mx program, modb.mx shows an alternative approach to the problem. Rather than 
pre-compute the moderator* predictor variable and use the multiple regression model, we use 
special features in Mx to place the observed moderator values directly on a path in the diagram. 
This time we are still regressing the outcome variable on the predictor and moderator main 
effects, but a second copy of the moderator variable is used as a path in the structural equation 
model. During model-fitting, Mx is recomputing the predicted covariance for each subject, so 
tliat every subject has a different model based on their moderator variable value. Again we fit the 
model with and without the moderator effect, and we recapture the same difference in 
log-likelihood between the two runs with 1 d.f. (i.e., 3698.475 - 3577. 86S = 120.6). 

^rhis method takes quite a bit longer to run than the first approach, so what is its value? 
Consideration of the further possibilities of this method gives the answer. We could very easily 
revise the script to employ some arbitrarily complex function of the moderator variable. Perhaps 
we expect some exponential model? Or a sinusoidal one? Whatever theoretical expectations we 
have, they can be modelled directly and uniquely for each subject. 



Moda.mx Program 

moda.mx: 

Title Moderated regression with raw data 

Data Ninput=4 Nobs=500 Ngroups=l 

Rectangular file=moda.raw ! Data file contains 4 variables: 

Labels Predictor Moderator Predictor*Moderator Outcome 

Matrices 

R symm 3 3 Free !Cov matrix of Pred, Mod, Pred*Mod 
L full 1 3 Free ! Regression of Outcome on above 3 variables 
D diag 1 1 Free ! Error of outcome 

M full 1 4 Free ! Matrix for estimated means of all 4 variables 
Means M; [Predicted means 
Cov R I R*L' , 

L*R I L*R*L'+D ; IPredicted covariance structure 
Start .5 all ! Starting values 
Start 3R1 1R22R33 

Option multiple-fit !To allow fitting of submodel after the end statement 
End 

Drop y !Fix at zero the 9th parameter, regression of outcome on pred*mod 
End 

Modb.mx Program 

modb.mx: 

Title Moderated regression with raw data 

#define nx 2 ! Number of X variables 

#define ny 1 ! Number of Y variables 

Data Ninput=4 Nobservations=5(X) Ngroups=l 

Labels Predictor ModeratorMainEffect Moderator Outcome 

Rectangular file=modb.raw 

Definition Moderator / 

Matrices 

R symm nx nx Free !Cov matrix of Predictor and ModeratorMainEffect (X vars) 
J full ny nx Free IRegression of Outcome on X variables 
K full ny nx !To contain individual moderator values for each subject 
L full ny nx IFree parameter for magnitude of moderating effect 
D diag ny ny Free lError of Outcome (Y) 
M full 1 3 Free !For estimated means of X & Y 
Means M; ! Matrix expression for means 
CovRIR*(J+L.Ky_ 
(J+L.K)*R I (J+L.K)*R*(J+L.K)'+D ; IMatrix expression for covariance 
Icodes 1 2 4 99 !To indicate that variables K 2 & 4 are being analyzed 
Specify K -1 0 ! Individual moderator values will go in K 
Specify L 100 0 IFree parameter for size of moderator effect 
Start .5 all 
Start .0 L 1 1 

Start 1 R 1 1 R 2 2 IStarting values 

Option Multiple-fit !To allow fitting of submodel 

End 

drop -1 100 ! Refit model with moderator effect fixed at zero 
End 

moda.raw: <insert raw data file here from SAS program via ftp address> 
modb.raw: <insert raw data file here from SAS program via ftp address> 
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OUTPUT OF MODA.MX 

(abridged) 



** Mx startup successful ** 

** MX-Windows version 1.30** 

Summary of VL t'iie data for group 1 

Code 1.0000 2.0000 3.0000 4.0000 

Number 500.0000 500.0000 500.0000 500.0000 

Mean 0.0314 -0.0189 0.5145 0.1258 

Variance 1.0162 1.0074 1.3608 1.0690 

MX PARAMETER ESTIMATES 

GROUP NUMBER: 1 

TITLE MODERATED REGRESSION WITH RAW DATA 
MATRIX D 

This is a DIAGONAL matrix of order 1 by 1 
1 

1 0.3392 
MATRIX L 

This is a FULL matrix of order 1 by 3 

1 2 3 
1 0.4972 0.4281 0.2609 

MATRIX M 

This IS a FULL matrix of order 1 by 4 

12 3 4 
1 0.0314-0.0189 0.5145 0.1258 

MATRIX R 

This is a SYMMETRIC matrix of order 3 by 3 
1 2 3 

1 1.0162 

2 0.5151 1.0074 

3 -0.0420 -0.0318 1.3608 

Your model has 14 estimated parameters and 2{)()() Observed statistics 
-2 times lug-likelihuud of data = 5150.144 
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Multiple Fit option in effect 

The foUowiag MX script lines have been read: 

DROP 9 
END 



Summary of VL file data for group 1 

Code 1.0000 2.0000 3.0000 4.0000 

Number 500.0000 500.0000 500.0000 500.0000 

Mean 0.0314 -0.0189 0.5145 0.1258 

Variance 1.0162 1.0074 1.3608 1.0690 



MX PARAMETER ESTIMATES 
GROUP NUMBER: 1 

TITLE MODERATED REGRESSION WITH RAW DATA 
MATRIX D 

This is a DIAGONAL matrix of order 1 by 1 
1 

1 0.4317 
M/xTRIX L 

This is a FULL matrix of order 1 by 3 

I 2 3 

I 0.4883 0.4245 0.0000 

MATRIX M 

This is a FULL matrix of order 1 by 4 

12 3 4 
1 0.0314 -0.0189 0.5145 0.1258 

MATRIX R 

This is a SYMMETRIC matrix of order 3 by 3 
1 2 3 

1 1.0162 

2 0.5151 1.0074 

3 -0.0420 -0.0318 1.3608 

Your model has 13 estimated parameters and 2()()() Observed statistics 
•2 times lug-likelihuod uf data = 5270.752 
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OUTPUT OF MODB.MX 



** Mx startup successful ** 
**MX -Windows version 1.30** 

Summary of VL file data for group' 1 

Code -1.0000 1.0000 2.0000 4.0000 
Number 500.0000 500.0000 500.0000 500.0000 
Mean -0.0189 0.0314 -0.0189 0.1258 
Variance 1.0074 1.0162 1.0074 1.0690 

MX PARAMETER ESTIMATES 

GROUP NUMBER: 1 

TITLE MODERATED REGRESSION WITH RAW DATA 



MATRIX D 

This is a DIAGONAL matrix of order 1 by 1 
1 

I 0.3392 
MATRIX J 

This is a FULL matrix of order 1 by 2 

1 2 
1 0.4972 0.4363 

mat:<ix k 

This is a FULL matrix of order 1 by 2 

1 2 
1 0.7378 0.0000 

MATRIX L 

This is a FULL matrix of order 1 by 2 

1 2 
1 0.2609 0.0000 

MATRIX M 

This is a FULL matrix of order 1 by 3 

1 2 3 

1 0.0314-0.01X9-0.00X6 

MATRIX R 

This is a SYMMETRIC matrix of order 2 by 2 
1 2 

1 1.0162 

2 0.5151 1.0074 

Your model has 10 estimated parameters and 2()()() Observed statistics 
•2 times log-likelihood of data = 3577.868 
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Multiple fit option in effect 

The following MX script lines have been read: 



DROP-1 100 
END 



Summary of VL file data for group 1 

Code -1.0000 1.0000 2.0000 4.0000 

Number 500.0000 500.0000 5()i).0()0() 5()().()()()0 

Mean -0.0189 0.0314 -0 0189 0.1258 

Variance 1.0074 1.0162 1.0074 1.0690 



MX PARAMETER ESTIMATES 
GROUP NUMBER: ) 

TITLE MODERATED REGRESSION WITH RAW DATA 



MATRIX D 

This is a DIAGONAL matrix of order I by 1 
I 

1 0.4317 



MATRIX J 

This is a FULL matrix of order 1 by 2 

1 2 
1 0.4883 0,4245 

MATRIX K 

This is a FULL matrix of order I by 2 

I 2 
1 0.7378 0.0000 



MATRIX L 

This is a FULL matrix of order 1 by 2 

1 2 
1 0.0000 0.0000 



MATRIX M 

This is a FULL matrix of order 1 by 3 

1 2 3 

I 0.0314 -0.0189 0.1258 



MATRIX R 

This is a SYMMETRIC matrix of order 2 by 2 
1 2 

1 1.0162 

2 0.5151 1.0074 

Your model has 9 estimated parameters and 2000 Observed statistics 
•2 times log-likelihood of data = 3698.475 
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AMOS PROGRAM EXAMPLE 



! nonlin.amd file contents 

! Nonlinear example data from Kenny and Judd. 1984 

SSample size = 500 
$input variables 
xl 

x2 
xlxl 
x2x2 
xlx2 

y 

SCovariances 

1.150 

.617 .981 
-.068 -.025 2.708 

.075 .159 .729 1.717 

.063 .065 1.459 1.142 1.484 

.256 .166-1.017 -.340 -.6 10. 763 



Note: The nonlin.amw file used in Amos to compute the coefficients is generated from the model 
specification diagram on the next page. The unstandardized coefficients are then printed 
for each path (unstandardized estimates) on the diagram (see subsequent page) 
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SAS/PROC CALIS EXAMPLE 



Interactive models don't often converge and are often difficult to specify. A two-stage fit seems to work 
best. A short SAS program is provided which will generate the interactive model, given the LINEQS and 
STD statements of the original measurement model with contraints on the covariances rather than 
generating nuisance variables. These covariance models run quicker and yield identical results to the 
model proposed by Kenny & Judd (1984). 

I have provided two SAS programs below: (1) the first fits the Kenny & Judd model using LISREL and 
the appropriate constraints: v2) the second is a small program which the user can employ by putting the 
LineEQS statements on one card section and the STD statements on another for the measurement model 
of interest. A punch file is then created which is the Proc CALIS program you need with all of the 
necessary constraints under a restatement of Kenny & Judd's model with only covariances. 

*This is the interactive example on page 208 of Loehlin's Latent Variable 
Modeling book; 

datacorrs (type=cov);input xl x2zl z2xlzl xlz2x2zl x2z2 y _name_ $ 58-61; 



_type_='COV';cards; 

2.39^ XI 

1.25^ 1.542 X2 

.445.202 2.097 Zl 

.231 .116 1.141 1.370 Z2 

-.367 -.070 -.148 -.133 5.669 .... XlZl 
-.301 -.041 -.130 -.117 2.868 3.076 . . . X1Z2 
-.081 -.054 .038 .037 2.989 1.346 3.411 . . X2Z1 



-.047 -.045 .039 -.043 1.341 1.392 1.719 1.960 . X2Z2 
-.368 -.179 .402 .282 2.556 1.579 1.623 .971 2.174 

Ydata size (type=cov):input xl x2 zl z2xlzl xlz2x2zl x2z2 y _name_ $58-61: 
cards; 

0 0 0 0 0 0 0 0 0 MEAN 

2.395 1 .542 2.097 1 .370 5.669 3.076 3.4 1 1 1 .960 2. 1 74 STD 

500 500 500 500 500 500 500 500 500 N 

data corrs (type=cov);set corrs size; 
proc print, 

proc calis cov data=corrs method=gls maxiter=300 tech=quanew 
edf=499; 
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Lineqsxl = f 1 + el ,x2 =g (.646) fl + e2 ,zl = f2 + e3 ,z2 =h (.685) f2 + e4 ,xlzl = flf2 +ele3 
,xlz2 =h f lf2 +ele4 ,x2zl =g flf2 +e2e3 ,x2z2 =gh f lf2 +e2e4 ,y=c f 1 + d f2 + e f lf2 + e5;stdel =erl 
,e2 =er2 ,e3 =er3 ,e4 =er4 ,e5=er5,ele3 =erler3 ,ele4 =erler4 ,e2e3 =er2er3 ,e2e4 =er2er4 ,f 1 =vf 1 ,f2 
=vf2 ,f lf2 =vflf2 ;covele3 ele4 =cl ,ele3 e2e3 =c2 ,ele4 e2e4 =c3 ,e2e3 e2e4 =c4 ,fl f2 =cf lf2 
;erler3 =erl *er3 +vfl *er3 +vf2 *erl ;erler4 =erl *er4 +vf 1 *er4 +h *h *vf2 *erl ;er2er3 =er2 *er3 
+g *g *vf 1 *er3 +vf2 *er2 ;er2er4 =er2 *er4 +g *g *vf 1 *er4 +h *h *vf2 *er2 ;vf lf2 =vf 1 *vf2 +cf lf2 
*cf lf2 ;cl =vf2 *erl *h ;c2 =vfl *er3 *g ;c3 =vfl *er4 *g ;c4 =vf2 *er2 *h ;gh =g *h ; 

♦Note: all variable names and parameters must be limited to 4 characters, because this program generates 
product variables and parameters which are products/functions of other variables. 8 character 
namesparameters will result in generation of illegal SAS names/values; 

data try;infile cards eof=last; 

array name{ 11}$ xvarl-xvarl l;array load{ 11}$ xloadl-xloadl l;array 

eld{ 1 1 } $ eloadl-eloadl l;array enmt' 11}$ enamel -enamel l;retain xvarl-xvarl 1 xloadl-xloadl 1 
eloadl-eloadl ] enamel -enamel 1 firsty lastx bottom factln fact2n;input vamame $ symbol $ fload $ 
fname $ symbol $ eload $ ename $; 

if fname=lag(fname) then cut=4);else cut=l;i-t-l; 

if cut=l then do;factln=lag(fname);fact2n=fname;lastx=i-l;firsty=i;end; 

name{i}=vamame; 

load{i}=fload; 

eld{i}-eload; 

enm{i}=translate(ename; ',','); 
ifeld{i}='rthen eld{i}="; 

if load{i}= r then load{i}= ';bottom=i;return;last:nvar=i;output; 
cards; 

vara = 11 fl + esdl eUvarb = 12 f 1 esd2 o2,varc = 13 fl -t- esd3 e3,vard 14 f 1 + esd4 e4,varg = ml f2 
+ esdS e5,varh = m2 f2 esd6 e6,vari = m;. f2 + esd7 e7,data try2; 
infile cards eof=last;array env{ 11}$ evar'-evarl 1; 
retain enamel -ename 1 1 evarl-evarll nvar2; 

input ename $ equals $ evar $ ;l-t-l;env(i}=translate(evar, ',','); if env{i| = T then env{i}="; 
return; 

last:nvar2=i;output; 
cards; 

el = l,e2 =: l,e3 = Ke4= Ue5 = l,e6 = Ke'^ = Kfl = Kf2 = Kdata try;set try: 
drop cut vamame symbol fload fname eload ename;da^ i try2;set try2; 
drop 1 ename equals evar;data try;merge ty try2;data try;set try; 
file punch; 

'iXmy name{ 11}$ xvarl-xvarl 1; 
an ay load{ 11}$ xloadl-xloadl 1; 
array eld { 1 1 } $ evar 1 -evar 1 1 ; 
array eld2 {11}$ eload 1 -eload 1 1 ; 
array enm{ 11}$ enamel -enamel 1; 
do i=l to lastx; do j=firsty to bottom; 
intvar=compress(name { i } llname { j } ); 
put intvar =' name{i} name{j} V; 
end;end; 

put Proc Calis cov method=gls;'; 

put *Lineqs';do i=l to nvar; 

put nameii} =' load{i} @; 

if i<firsty then put factl n @;else put fact2n @; 
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put V@; 

ifeld2{i}>"then puteld2{i} " @; 

putenm{i} \';end;do i=l to lastx; 

do j=firsty to bottom; 

int\'ar=coinpress(name{i }ttname{ j } ); 

intload=compress(load{ i }llload{ j } ); 

intfaame^scompressCf act 1 nllf act2n); 

intenamc=comprcss(cnm{ i } llenm { j } ); 

put intvar intload Intfname intename @; 

if i=lastx and j=bottom then put V;else put 

end;end;put 'std';do i=l to bottom;if eld{i}=' ' then eld{i}= T; 

ifeld2{i}="then eld2{i}='r;put enm{i} eld{i} \';end; 

do i=l to lastx; 

do j=firsty to bottom; 

intename=compress(enm{ i } tlenm{ j } ); 

if eld { i } >• 1 • then texta=eld { i } ;else texta=eld2 { i } ; 

if eld{j}>'r thtn textb=eld{j};else textb=eld2{j}; 

inteload=compress(textalttextb); 

put intename inteload \'; 

end;end; 

put fact In W fact In ';; 
put fact2n '=v* fact2n 
put intfname =v' intfname 
put 'cov';ind=0;do i=l to lastx; 
do j=firsty to bottom; 
intename=compress(enm{ i }ttenm{ j } ): 
do k= i to lastx: 
do l=j to bottom; 

intenam2=compress(enm{ k } llenm { 1 } ); 
if (i=k or j=l) and not(i=k and j=l) then do; 

ind=ind+l; 

cname=c1tind; 

cname=compress(cname); 

put intename ' ' intenam2 cname \': 

end; end; end; end; end; put factln ' ' fact2n '=c' intfname V; 

"^here are the covariance constraints; 

do i=l to lastx; 

do j=firsty to bottom; 

if eld { i } >' r then texta=eld { i } ;else texta^eld2 { i } ; 
if eld{j}>'r then textb=eld{j};else textb=eld2{ j}; 
intename=compress(textalltextb); 
put intename '= @; 
ifeld{i}yrthen puteld{i} @; 
if eld2 { i }>• r then put eld2 { i }'* eld2 { i } (S) ; 
put @; 

if eld { j } >• r then put eld { j } @ ; 

if eld2{ j }>' 1 ' then put eld2{ j } eld2{ j } @; 

put V®; 

if load { i } >" then put load { i } '* ' load { i } @ ; 
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put Vfactln***®; 
ifeld{j}>'rthen puteld{j} V@; 
ifeld2{j}>'r then puteld2{j} •*'eld2{j} V@; 
ifload{j}>"then putload{j} •*'load{j} '*'@; 
put Vfact2n '*'(a); 
ifeld{i}>'rthen puteld{i} V; 
ifeld2{i}>'rthen puteld2{i} •*'eld2{i} V; 
end;end; 

intfname=compress(factlnllfact2n); 

put V intfname =v' factln '*v* fact2n '+c' intfname intfname V;ind=(); 

^constraints for the covariances; 

do i=l to lastx; 

do j=firsty to bottom;if eld{i}>'r then intename=compress(eld{i}lleld{j});if eld2{i}>'r then 
intename=compress(eld2{ i }lleld{ j }); 

do k= i to lastx; 

do l=j to bottom; 

if eld{k}>*r then intename=compress(eld{k|lleld{l}); 
if eld2{k}>'r then intename=compress(eld2{k}lleld{l}); 
if (i=k or j=l) and not(i=k and j=l) then 
do; 

ind=ind+l; 
cname=c'llind; 
cname=compress(cname); 
putciiame ='@; 
if i=k then do; 

put V fact2n '*'(§); 

ifeld{i}>'rthen put eld{i} @; 

if eld2 { i } >' r then put eld2 { i } eld2 ( i } @ ; 

put load{l} V; 

end; 
if j=l then do; 

put V factln'*'@; 

ifeld{j}>'rthen put eld{j}@; 

if eld2{j}>'r then put eld2{j} eld2{j} @; 

put load{k} V; 

end; 

end; end; end; end; end; 
do i=l to lastx; 
do j=f irsty io bottom; 

if load{i}>' ' and load{j}>' 'then do; 

intfname=load { i } llload { j } ; 

intfname=compress(intfname); 

put intfname =' load{i} load{j} V; 

end; end; end; 
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